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ABSTRACT Paramyxoviruses, a family of RNA enveloped viruses that includes human parainfluenza virus type 3 (HPIV3), cause 
the majority of childhood croup, bronchiolitis, and pneumonia worldwide. Infection starts with host cell receptor binding and 
fusion of the viral envelope with the cell membrane at the cell surface. The fusion process requires interaction of the two viral 
surface glycoproteins, the hemagglutinin-neuraminidase (HN) and the fusion protein (F). We have previously shown that vi- 
ruses with an HN/F pair that is highly fusogenic in monolayers of immortalized cells due to mutations in HN's secondary sialic 
acid binding site are growth impaired in differentiated human airway epithelium (HAE) cultures and in vivo. Here we have 
shown that adaptation of HPIV3 to growth in the lung is determined by specific features of HN and F that are different from 
those required for growth in cultured immortalized cells. An HPIV3 virus bearing a mutated HN (H552Q), which is fit and fuso- 
genic in immortalized cells but unfit for growth in the lung, evolved into a less-fusogenic but viable virus in differentiated hu- 
man airway epithelium. Stepwise evolution led to a progressive decrease in efficiency of fusion activation by the HN/F pair, with 
a mutation in F first decreasing the activation of F by HN and a mutation in HN's secondary sialic acid binding site decreasing 
fusion activation further and producing a stable virus. Adaptation of HPIV3 to successful growth in HAE is determined by spe- 
cific features of HN and F that lead to a less easily activated fusion mechanism. 

IMPORTANCE Human parainfluenza viruses (HPIVs) are paramyxoviruses that cause the majority of childhood cases of croup, 
bronchiolitis, and pneumonia worldwide, but there are currently no vaccines or antivirals available for treatment. Enveloped 
viruses must fuse their membrane with the target cell membrane in order to initiate infection. Parainfluenza virus fusion pro- 
ceeds via a multistep reaction orchestrated by the two glycoproteins that make up its fusion machine. The receptor-binding 
hemagglutinin-neuraminidase (HN), upon receptor engagement, activates the fusion protein (F) to penetrate the target cell and 
mediate viral entry. In this study, we show that the precise balance of fusion activation properties of these two glycoproteins dur- 
ing entry is key for infection. In clinically relevant tissues, viruses evolve to acquire a set of fusion features that provide key clues 
about requirements for infection in human beings. 
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Human parainfluenza virus 3 (HPIV3) is a major cause of pe- 
diatric croup, bronchiolitis, and pneumonia worldwide, and 
in the United States it accounts for a significant share of pediatric 
respiratory hospitalizations, as well as a rising number of infec- 
tions in immunocompromised patients. There are no vaccines or 
treatments against it or any of the other HPIVs (1-4). HPIVs are 
enveloped paramyxoviruses, carrying two surface glycoproteins, 
the hemagglutinin-neuraminidase (HN) and the fusion (F) pro- 
tein. After binding to the host receptor, sialic acid, HN and F 
interact to mediate fusion prior to delivery of the viral genome 
into the host cell. The correct timing of F activation by HN is 
essential for infection: to be effective, triggering of F to its active 
state must occur only when F is positioned near the target cell 
membrane (5). 

The HN protein, which is present as a tetramer composed of 
two dimers on the surface of the virus and the infected cell, per- 



forms three important functions during the HPIV3 viral life cycle. 
In addition to receptor binding, the initial step of entry, and F 
activation, which promotes the steps of entry after receptor bind- 
ing, HN has the third function of receptor cleavage, performed by 
its neuraminidase activity to promote the budding of new viruses 
from an infected cell. The communication between the two mol- 
ecules that comprise the paramyxovirus fusion machine is biolog- 
ically relevant, with impact on infection in the natural host. We 
have studied HPIV3 virus variants bearing mutations at one of 
HN's sialic acid receptor binding sites, either in the primary re- 
ceptor binding site identified in the existing crystal structure or in 
a secondary sialic acid binding site (site II), present at the dimer 
interface and required for activation of F (6-8). Alterations in site 
II at the dimer interface of HPIV3 HN influence growth, implicat- 
ing HN's site II in pathogenesis (9). The virus HPIV3 HN H552Q, 
bearing a mutation at HN's site II, binds its receptor more avidly, 
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activates F more efficiently, fuses in monolayer culture more ex- 
tensively (6, 7), and physically interacts with F more strongly (10). 
However, this virus is unfit in vivo (cotton rat model) and in hu- 
man airway epithelium (HAE), a pseudostratified, highly differ- 
entiated model of unpassaged primary cells that closely resembles 
the epithelial tissue of the human respiratory tract (9). This sug- 
gests that an increase in F activation efficiency is a disadvantage in 
the natural host. The virus with enhanced F triggering produces 
inactive viral particles in HAE, in which the F protein may have 
been triggered before contacting target cells (9). These HN-F re- 
lationships are physiologically relevant (9): too-active HN-F in- 
teraction may lead to premature, and therefore ineffective, activa- 
tion of the fusion process (5). 

To obtain experimental support or proof for this notion, we 
turned to the virus and the process of natural selection: in HAE, 
which includes a variety of cell types, cilia, and mucus, as in the 
human lung (11), the site II-mutated variant HPIV3 HN H552Q 
virus adapted for survival by evolving to decrease the ease of trig- 
gering of F and the extent of fusion. We have shown here that 
infection in lung tissue depends upon a specific balance of fusion 
activation, determined by both the efficiency of HN at triggering F 
and the stability of the pretriggered state of F. Viruses with highly 
efficient F-triggering capacity — whether HN triggers F too effi- 
ciently or F is too readily triggered — were not favored in the selec- 
tion process, likely due to premature triggering of F. The HPIV3 
HN H552Q variant virus under the selective pressure of HAE ac- 
quired sequential mutations in the F and HN proteins that low- 
ered HN-receptor interaction and lowered F activation, ultimately 
promoting growth in HAE. The first mutation emerged in the F 
protein, followed by a mutation in HN's secondary sialic acid 
binding site. No reversion to the laboratory wild-type (wt) virus 
was identified. The new HPIV3 variant that emerged during 
growth of HPIV3 HN H552Q in HAE carries an HN-F fusion 
apparatus with functional properties that are suited to growth in 
lung tissue. 

RESULTS 

Growth of HPIV3 HN H552Q in HAE is ineffective but leads to 
emergence of plaque size variant viruses. The variant HPIV3 HN 
H552Q infected HAE cells, but by 3 days it did not produce infec- 
tious progeny virions, yielding only noninfectious viral particles 
(9, 12) . To determine whether the selective pressure of HAE would 
yield a virus suited to growth in HAE, HPIV HN H552Q was 
grown in HAE for 12 days. Cultures were infected with 4,000 PFU 
of either HPIV3 wt or HN H552Q, fluids containing released vi- 
ruses were collected each day postinfection, and titers were deter- 
mined in cell monolayers (Fig. 1). 

Figure 1A outlines the experiment. Infection with HPIV3 wt 
(blue line) produced infectious titers throughout the experiment, 
with plaque size remaining unchanged after 12 days. In contrast, 
infection with HPIV3 HN H552Q (bright red) yielded zero viral 
titer up to 3 days after infection. At 5 days after infection, a mixed 
population (brown) of medium-plaque (diameter, 0.39 mm) and 
small-plaque (diameter, 0.12 mm) viruses was retrieved in repli- 
cate experiments, unlike either the large plaques (diameter, 
1.08 mm) of the highly fusogenic parental HPIV3 HN H552Q or 
plaques of HPIV3 wt (diameter, 0.61 mm), as shown in Fig. IB. By 
7 to 9 days after infection, the population contained only small- 
plaque viruses, at titers greater than those of HPIV3 wt (data not 
shown), and medium-plaque virus did not reappear. This pattern 



of emergence of variant viruses consistently reoccurred in each of 
the seven times that this experimental protocol was carried out. 

To assess the growth in HAE of the medium- and small-plaque 
viruses relative to that of parental HPIV3 HN H552Q, samples 
were plaque purified and used to infect HAE. Their growth was 
compared to those of the parent and wild-type viruses (Fig. 1C). 
As expected, up to 3 days postinfection, HPIV3 HN H552Q 
yielded zero titer, while the HPIV3 wt produced 10 5 PFU/ml. In- 
fection with the medium-plaque virus yielded up to 10 4 PFU/ml 
2 days postinfection, followed by emergence of a mixed popula- 
tion of medium and small plaques, with only small plaques by 
7 days postinfection (data not shown). Infection with the small- 
plaque virus yielded growth up to 3 logs higher than that of HPIV3 
wt within the first 3 days of HAE infection, peaking at 10 7 PFU/ml. 
The HN H552Q variant appears to have undergone sequential 
alterations that overcame its growth defect in HAE. The associ- 
ated, progressive decrease in plaque size (Fig. IB) suggests selec- 
tion in HAE for viruses with a lower cell fusion capacity. 

The high infectious titers attained in HAE by the small-plaque 
virus relative to those of the wt virus, as well as HPIV3 HN H552Q, 
might be accounted for by an increase in the number of released 
virions (total particles). However, quantification by quantitative 
real-time PCR (qRT-PCR) of the amount of genome retrieved 
from HAE 3 days after infection for the small-plaque virus was 
similar to that for wt virus of HPIV3 HN H552Q (Fig. ID). This 
similarity permitted calculation of the amounts of infectious viri- 
ons. On average, 60% of the small-plaque virus population was 
infectious, as opposed to 15% of the HPIV3 wt population and 0% 
for the HPIV3 HN H552Q variant (Fig. IE). Thus, while all three 
viruses produce particles, the ratio of infectious particles to total 
particles is higher in the stable small-plaque virus that emerges 
from the HPIV3 HN H552Q parent in HAE. 

Viability of emerging variants is associated with mutations 
in HN and F envelope glycoproteins. The RNA of the viruses 
emerging from six different infections of HAE with HPIV3 wt, HN 
H552Q, the medium-plaque viruses, and the small-plaque viruses 
was extracted, and the HN and F genes were sequenced (Table 1). 
Viruses retrieved from HAE infected with HPIV3 wt retained wt 
HN and F sequences. Noninfectious HPIV3 HN H552Q collected 
1 day postinfection carried HN with the H552Q mutation and 
wild-type F. The medium-plaque virus (the first to appear) carried 
HN H552Q and F G396D. The small-plaque virus carried the HN 
H552Q and F G396D mutations, as well as a Q559R mutation in 
HN. Figure 2A shows the HPIV3 HN dimer structure, with 
H552Q and Q559 marked; Fig. 2B shows the homology model of 
the HPIV3 F protein prefusion structure, and the HPIV3 F post- 
fusion structure, with G396 marked. The Q559R HN mutation is 
located in the secondary sialic acid binding site that is present at 
the dimer interface and has been implicated in activation of F (7, 
13). 

Evolution towards slower HN-F triggering and fusion. The 

failure of the highly fusogenic HPIV3 HN H552Q virus to produce 
active virions in HAE infection could be due to the premature 
triggering of the F protein in the progeny virions, leading to inac- 
tivation of the virions (9). In this case, positive selection in HAE 
should favor mutations that overcome this tendency to prema- 
turely trigger F. A quantitative fusion assay was used to evaluate 
the fusion promotion capacity of each HN-F pair (7, 13, 32). 293T 
cells coexpressing F wt or F G396D with HN H552Q or HN 
H552Q,Q559R were overlaid with receptor-bearing cells at 37°C 
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FIG 1 HAE infection with HPIV3 HN H552Q: emergence of variants. (A) Timeline of viral growth after infection of HAE with wt HPIV3 (blue) and HPIV3 HN 
H552Q (red). (B) Plaque diameter in CV1 cells 16 to 18 h after infection with the indicated virus retrieved from HAE (n = 20 per virus from 3 experiments). (C) 
Titers of viruses grown in HAE: wt HPIV3 (blue), HPIV3 HN H552Q (red), or HPIV3 "small plaques" (brown). Data points are means (±SD) of triplicate 
measurements and are representative of at least 5 experiments. (D) Quantitation of viral genome released 3 days post-HAE infection (determined by qRT-PCR). 
Data points are means (±SD) of triplicate experiments. (E) Percentage of infectious virions (out of total particles) retrieved from infection of HAE. 



for 3 h, and fusion was measured. As observed with plaque sizes 
(Fig. IB), fusion is reduced if F G396D (instead of wt F) is paired 
with HN H552Q. Fusion decreased by 32% when wt F is replaced 



with F G396D, indicating that the mutation G396D produced less 
fusion in the presence of the same H552Q HN (Fig. 3A). With the 
second HN mutation (H552Q Q559R), the fusion capacity de- 
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TABLE 1 Mutations acquired in HN and F glycoproteins of HPIV3 
variants 



HPIV3 variant 


HN mutation(s) 


F mutation 


HN H552Q 


H552Q 


None 


Medium-plaque virus 


H552Q 


G396D 


Small-plaque virus 


H552Q, Q559R 


G396D 



creases more, indicating that the Q559R mutation in HN also 
decreased the efficiency of the HN-F fusion pair. 

Specific HN-F-triggering efficiency was assessed using a 
fusion-triggering assay to determine the temperature of F activa- 
tion for each HN-F pair. Cells coexpressing the indicated HN and 




Pre-fusion F Post-fusion F 

FIG 2 Molecular structure of wild-type and mutant HN and F. (A) HN dimer 
structure (PDB IV2I), with the monomers colored gray and cyan. The loca- 
tions of the H552Q and Q559R mutations at the dimer interface are high- 
lighted by showing residues H552 and Q559 as sticks and marked with aster- 
isks. The structure was subjected to 20 ns of MD simulation. (B) Crystal 
structure of the postfusion HPIV3 F protein (right) and a homology model of 
the prefusion protein, based on PIV5 F (left). The three chains of F are colored 
green, cyan, or magenta. The G396 residues are represented as yellow spheres. 



F pair (Fig. 3B) were allowed to bind erythrocytes (RBCs) at 4°C. 
The cells were then washed to remove unbound RBCs and trans- 
ferred to fresh medium at 4, 15, 22, 30, or 37°C for 1 h, to allow F 
activation. The amount of target RBCs that (i) were released into 
the medium by HN's innate neuraminidase activity (dotted line), 
(ii) were bound but had not fused (dashed lines), or (iii) had 
undergone fusion (solid lines) was determined. The temperature 
of F activation (defined as the temperature for 50% RBC fusion) 
increased from 12°C for HN H552Q with F wt to 29°C for HN 
H552Q with F G396D (medium plaques), indicating that the F 
mutant requires more heat energy to be activated than does F wt in 
the presence of the same HN (Fig. 3B). The virus yielding stable 
and infectious small plaques, bearing HN H552Q,Q559R and F 
G396D, had the highest temperature requirement (>37°C) for F 
activation, achieving only 20% fusion at 37°C. The final extent of 
fusion accomplished by this HN/F pair is also affected by the sig- 
nificant amount of release of RBCs from the expressing cell sur- 
faces; this issue will be further addressed below (see Fig. 5C). 

The mutation G396D confers slower activation kinetics of 
the F protein. The effect of the G396D mutation in F could result 
from less-efficient triggering or from less-effective fusion promo- 
tion after triggering. To distinguish between these two possibili- 
ties, heat activation of wt F and mutant F G396D in the absence of 
HN was carried out. Hemagglutinin (HA) can be used to bind to 
receptor-bearing cells that are targets for HPIV3 F-mediated fu- 
sion; since influenza virus HA and HPIV3 F do not interact on the 
cell surface and HA does not trigger F to mediate fusion (10), HA 
serves simply to cause the cells to adhere. Cells were transfected 
with influenza virus HA with either wt F or mutant F G396D, the 
cells were incubated at 45°C to activate F ( 14, 15), and fusion was 
measured over time (Fig. 4). At 30 min, F G396D showed approx- 
imately 18% less fusion than F wt. By 60 min there was no signif- 
icant difference in the fusion between F wt and F G396D, indicat- 
ing that F G396D was less efficient at initial activation but once 
activated was fully fusion competent. 

Mutations acquired in HN's secondary sialic acid binding 
site ("site II") lower avidity. To dissect the role of the HN Q559R 
mutation in the secondary sialic acid binding site, HNs bearing 
both H552Q and Q559R were compared with each singly mutated 
HN and wt HN for receptor avidity, neuraminidase, and 
F-triggering properties. To measure receptor binding avidity, cells 
transiently expressing either of the four HNs were pretreated with 
neuraminidase to remove receptors on the expressing cells' sur- 
faces and eliminate potential interactions between HN and recep- 
tor molecules on the HN-expressing cell surface (7). Aliquots of 
RBCs with graded receptor depletion were then used to quantify 
RBC binding to the HNs; HN molecules with higher avidity bind 
RBCs that have lower receptor density (higher depletion). As ex- 
pected, HN H552Q showed higher avidity than HN wt (50% bind- 
ing to RBCs treated with 54 mU of neuraminidase, compared to 
50% binding to RBCs treated with 44 mU neuraminidase) 
(Fig. 5A). For HN H552Q.Q559R, avidity is markedly reduced 
(50% binding at 19 mU neuraminidase). This decrease can be 
accounted for by the Q559R mutation, which singly (HN Q559R) 
showed 50% binding at 17 mU of neuraminidase. The strikingly 
reduced avidity conferred by Q559R points to the role of the sec- 
ondary sialic acid binding site in receptor avidity. 

The Q559R mutation also impacted HN's neuraminidase ac- 
tivity, known to reside in site I. Neuraminidase activity of the 
doubly mutated HN H552Q.Q559R was approximately 80% of 
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FIG 3 Cell fusion promoted by parental and mutant glycoproteins. (A) Cell-cell fusion mediated by the indicated HN-F pairs was quantitated after 3 h at 37°C 
using a beta-galactosidase complementation assay. The bars represent means ± SD of data from at least 3 experiments. (B) 293T cells coexpressing the indicated 
HN-F pair were allowed to bind RBCs at 4°C and then transferred to 4, 15, 22, 30, or 37°C for 1 h. Released, bound, and fused RBCs were quantified. The values 
are means (±SD) of results from three experiments. 



that of parental HN H552Q (see Table 2) . This is likely to be due to 
an allosteric effect on site I, but this remains to be explored mech- 
anistically. Since both the binding and neuraminidase activities of 
HN were affected by the Q559R mutation, we hypothesized that 
the duration of receptor engagement may be shorter for the small- 
plaque virus' HN than for parental HN H552Q. To address this 
question, the time course of each HN's binding and releasing sialic 
acid receptors was determined at 37°C. Cells expressing HN WT, 
HN H552Q, HN H552Q,Q559R, or HN Q559R were incubated 
with 1% RBC suspensions for 30 min at 4°C, rinsed to remove 
unbound RBCs, and transferred to 37°C. At the indicated time 
points (Fig. 5B), the released and bound RBCs were collected and 
quantified. Within 15 min of incubation at 37°C, approximately 
44% and 47% of RBCs were released from HN H552Q,Q559R and 
HN Q559R, respectively, compared to only 13% for HN wt and 
0.7% for HN H552Q. This faster receptor detachment from HN 
H552Q.Q559R and HN Q559R reflects the lower ratio of receptor 
avidity to neuraminidase activity for these variant HNs than for 
parental HN H552Q. 

Mutations acquired in HN's site II impair HN's F-triggering 
ability. To address the effect of Q559R on HN's F activation prop- 
erty, the fusion-triggering assay was carried out with HN wt, HN 
H552Q,Q559R, or HN Q559R coexpressed with F wt. Again the 




Time (minutes) 

FIG 4 Activation kinetics of F wt and F G396D in the absence of HN at 45°C. 
293T cells transiently expressing influenza virus HA and either F wt or F 
G396D were overlaid with RBCs at 45°C for 60 min. The percent RBCs fused 
(y axis) over time (x axis) was quantified. The values are means (±SD) of 
results from three experiments. 
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FIG 5 Effects of HN mutation on avidity, neuraminidase, receptor release, and F activation. (A) A panel of RBCs with different degrees of receptor depletion 
was used to quantify HAD on 293T cell monolayers expressing the indicated HNs. The values are means (±SD) of results from three experiments. (B) 293T cells 
transiently expressing the indicated HN were allowed to bind RBCs at 4°C for 30 min at pH 5 and then transferred to 37°C. Percent RBCs released at 0, 5, 10, 15, 
30, and 60 min were determined (y axis). Values are means (±SD) of results from triplicate wells. (C) 293T cells coexpressing the indicated HN-F pair were 
allowed to bind RBCs at 4°C and then transferred to 4, 15, 22, 30, or 37°C for 1 h. RBCs released, bound, and fused were quantified. Values are means ( ±SD) of 
results from three experiments. 



temperature of F activation was considered to be the temperature 
at which there is 50% RBC fusion. HN wt required a temperature 
of 18°C to achieve 50% binding or fusion (Fig. 5C). HN H552Q, 



the most efficient at F triggering, required only 12°C (Fig. 3B). 
Introduction of Q559R into the parental HN H552Q increased the 
activation temperature to 31°C, while singly mutated HN Q559R 



TABLE 2 Characteristics of variant HPIV3 HN and F proteins 



Protein 


Receptor avidity 
(mU for 50% 
RBC binding) 


Neuraminidase 
level 

(% of HN 
wt level) 


Release 

(% RBC released 
in 15 min at 
pH5) 


Temp of F activation (°C) 
F wt F G396D 


HN wt 


44 


100 


13 


18 


37 


HN H552Q 


54 


79.5 


0.7 


12 


29 


HN H552Q.Q559R 


19 


62.2 


44 


31 


>37 


HN Q559R 


17 


34.4 


47 


37 


>37 
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required a temperature of 37°C for 50% binding/fusion. The 
Q559R mutation confers a decrease in the efficiency of HN's 
F-triggering property. 

It is conceivable that the F activation temperatures of HN 
H552Q,Q559R and HN Q559R might have been obscured by the 
significant release of RBCs during the fusion-triggering assay. To 
address this issue, we removed the neuraminidase activity of these 
HNs in order to evaluate simply the inherent F-triggering ability 
of the HNs. The mutation D2 16R, which abolishes neuraminidase 
activity without affecting binding or F activation (13), was intro- 
duced into HN H552Q.Q559R and HN Q559R. HN carrying the 
D216R mutation is irreversibly receptor engaged. The HN 
D216R,H552Q,Q559R construct did not produce a functional 
protein on the cell surface, but HN D216R,Q559R was expressed 
at levels comparable to those for all other HNs (data not shown). 
Using HN D216R,Q559R in the presence of F wt, the fusion- 
triggering assay revealed that even without release of RBCs during 
the fusion assay, the temperature of activation for HN 
D216R.Q559R was still 28°C (Fig. 5C). This suggested that HN 
Q559R was inherently less efficient at activating F relative to HN 
wt and HN H552Q, a result not explained simply by release of 
receptor-bearing cells prior to F activation. 

Coevolution of HN and F in HAE identifies key characteris- 
tics required for HPIV3 growth in the lung. To achieve optimal 
growth characteristics in host tissue, the large-plaque HPIV3 HN 
H552Q virus underwent evolution to reduce the efficiency of its 
fusion machinery in several ways: by decreasing F's susceptibility 
to activation, by decreasing HN's receptor avidity, by decreasing 
HN's neuraminidase, and by decreasing HN's fusion-triggering 
capacity. The combination of characteristics led to diminished 
HN-F activation. Table 2 summarizes the individual characteris- 
tics of each protein in the series. The parent HPIV3 HN H552Q 
first acquired a mutation in F that impacted F's ability to be acti- 
vated, but this mutation alone was insufficient for generating a 
stable virus. The introduction of Q559R into HN's site II resulted 
in a decrease in avidity, neuraminidase, and HN's F-triggering 
functions, producing a virus that is viable in the lung tissue. 

DISCUSSION 

Infection by HPIV3 in lung is an ongoing positive selection pro- 
cess in which the strength of binding to host cells, fusion activa- 
tion, and receptor cleavage fluctuate until a stable balance is 
achieved. Our previous studies have shown that avid binding of 
HN to the host cell receptor, paired with efficient fusion activation 
by HN, was suited for infection and growth in monolayer culture 
but proved to be counterproductive in cotton rats and HAE (9). 
The studies reported here on the evolution of the virus in natural 
tissue, achieved through experiments that permit the HN/F fusion 
machinery to adapt in the natural host, explain the differences in 
outcomes in terms of specific HN-F activities over time. 

A key factor for viral infection in the host is an HN that does 
not prematurely activate F before viral entry. Several viral proper- 
ties contribute to this effect. In the lung, HPIV3 bearing an HN 
with lower avidity may be superior to one with a more avid HN 
because the resulting requirement for higher activation energy 
protects F from activation, delaying the triggering of F until the 
conditions are right for fusion. We have suggested (5, 9) that the 
correct timing of F activation is essential for viral entry and spread 
and that triggering should occur only when F is in proximity to the 
target cell membrane. HPIV3 HN H552Q produced noninfec- 



TABLE 3 Flow cytometric analysis of HN cell surface expression by 
FACS analysis with anti-HN monoclonal antibody 





No. of cell surface 


Protein 


HN molecules/cell" 


HN wt 


91,200 ± 24,617 


HN H552Q 


90,041 ± 19,838 


HN H552Q.Q559R 


93,263 ± 14,342 


HN Q559R 


96,878 ± 12,439 



" Data are representative of results from four experiments (averages ± SE). 



tious particles after primary infection of HAE (9), and we now 
have found that exposure to the selective pressure of growth in 
HAE favored the emergence of a new viral variant with slower 
fusion kinetics. It is possible that the generous availability of re- 
ceptors in the lung could be the factor in HAE that caused HPIV3 
HN H552Q viruses to be noninfectious after the primary infec- 
tion. HN H552Q, avid for receptor binding and efficient at F ac- 
tivation, would result in these viruses binding to sialic acid- 
containing moieties in the lung, leading to premature activation of 
F. The virus that emerged under the selective pressure in HAE 
(HPIV3 HN H552Q,Q559R F G396D) carries an HN with lower 
avidity for receptor binding and less-efficient activation of F. This 
variant grows even better than our lab-adapted wt HPIV3 in the 
HAE. These slow fusion machinery characteristics may also be 
shared with HPIV3 clinical isolates; an analysis is ongoing. Spe- 
cific fusion kinetics maybe a key factor in long-term virus survival 
and growth in vivo. 

The appearance and survival of the mutations in HN and F 
over time suggest several key characteristics of successful HN-F 
interactions. In host tissue, from a highly avid and fusogenic 
HPIV3 carrying HN H552Q, a new and unstable variant was pos- 
itively selected with a mutated fusion protein that was less fuso- 
genic than its predecessor. This unstable virus then acquired a 
second mutation in HN's secondary sialic acid binding site — 
Q559R — to form a stable and poorly fusogenic virus, which grew 
better in HAE than the HPIV3 wt virus that has been propagated 
in the laboratory. These findings suggest that the emergence of the 
less-fusogenic F protein mutant required changes in HN's site II to 
achieve a stable virus. Conversely, the mutations that were advan- 
tageous for growth in natural tissue led to poor growth in immor- 
talized cells; infection of CV1 cells with HPIV3 HN 
H552Q,Q559R/F G396D led to the acquisition of new mutations 
in the glycoproteins and a change in plaque size and morphology 
(data not shown). 

The H552 and Q559 residues are located at the HN dimer 
interface, near the second active site in the published X-ray crystal 
structure of HN (PDB 1V2I) (16) (Fig. 2A). The introduction of 
Q559R was the final step in the evolution from a virus unfit for 
growth in lung to a stable HPIV3 variant that is highly infectious in 
HAE. While the H552Q mutation confers higher avidity for recep- 
tors and increased F triggering relative to results for HN wt, the 
introduction of Q559R into the HN protein has the opposite ef- 
fect. Q559R decreases HN's avidity for receptors as well as its 
ability to activate the F protein. Note that the HN surface expres- 
sion levels do not differ between wt HN and mutants, including 
HN containing the Q559R mutation (Table 3). Using the pub- 
lished HN structure (above), the H552 residue was mutated to Q 
in silico. The mutation did not significantly alter the dimer inter- 
face or the local protein conformation. Based on this modeling 
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and molecular dynamics simulation, the 552Q mutation does not 
seem to disrupt and might in fact strengthen the interaction and 
participates in intermolecular hydrogen binding; this could ac- 
count for our recent biochemical data showing that H552Q alters 
dimerization of HN (10). 

The second mutation at this site, Q559R, might disrupt the 
dimer observed in the crystal structure. The arginines facing each 
other, in addition to not fitting into the tight space at the dimer 
interface, could create charge repulsion, possibly disrupting the 
dimer. However, generating this mutation in silico is not possible 
because the dimer interface of the currently available crystal struc- 
ture of HN (16) does not permit this mutation. Additional struc- 
tural analysis will explore the effects of this mutation on HN. 

The residue mutated in the F protein, glycine 396, is located on 
a long loop of F's domain II, which is surface exposed both in the 
postfusion crystal structure (PDB 1ZTM) (17) and in a homology 
model of the prefusion structure (model based on the parainflu- 
enza virus 5 [PIV5] F protein, PDB 2B9B) (18) (Fig. 2B). The 
G396D position is accessible and remains in the same structural 
location in the pre- and postfusion conformations of F. The mu- 
tation confers slower triggering kinetics on the protein, indicating 
that this domain is involved in the triggering mechanism. While in 
our study the mutation in F emerged to permit growth in lung 
tissue, several experiments using site-directed mutagenesis in 
other paramyxoviruses have shown that mutations in this domain 
of F can modulate F activation, most likely by stabilizing the pre- 
fusion state (19-21). Of note, this residue is the same location on 
F that was targeted by the effective prophylactic monoclonal an- 
tibody for respiratory syncytial virus, supporting the notion that 
this location is a key functional region that is preserved between 
paramyxovirus F proteins (22, 23). 

The progressive evolution from the parental HPIV3 HN 
H552Q to HPIV3 HN H552Q.Q559R/F G396D, from larger to 
smaller plaques and to decreased F triggering, highlights the dis- 
advantage of the parental virus' fusion machinery. The HN-F ma- 
chinery of HPIV3 HN H552Q, to adapt to the lung tissue, evolved 
to decrease its F activation capacity first, to grow in HAE, and a 
second mutation in HN's second sialic acid binding site, decreas- 
ing avidity and HN's F activation, conferred optimal growth in the 
airway cultures. The results support the role of HN's sialic acid 
binding site II in controlling F activation and underscore the no- 
tion that timely F activation is required for infection in host tissue. 
The HN-F pair of the HAE fit virus requires more energy to initi- 
ate fusion. In distinction from the features that confer fusogenicity 
in monolayer cultures, a virus carrying an HN with low receptor 
avidity, a high rate of release from the receptor, and decreased 
efficiency of activating the F protein and with reduced fusion is fit 
for growth in the natural host. 

The mechanisms underlying HN-F activity offer opportunities 
for subverting the viral fusion process. Apart from inhibiting viral 
entry, it maybe possible to force viral evolution towards develop- 
ing HN host receptor binding that is strong enough to prema- 
turely activate F and circumvent fusion. Recently, we identified 
compounds that stimulate HPIV3 HN to trigger F prematurely, 
before the virus contacts the target cell membrane, and thereby 
inactivate infectivity in HAE (5). Such compounds are of potential 
significance as antiparamyxovirus therapies and should be highly 
unlikely to elicit resistance mutations. Mutants arising under the 
selective pressure of premature-triggering compounds would 
have fusion machineries that are too impaired for viability and 



would thus be less fit and transmissible than the parent virus. 
Virions harboring mutations in their fusion machinery that would 
make them resistant to prematurely activating compounds would 
therefore be likely to be unresponsive in physiological situations, 
rendering them less fit to infect and cause disease. Characterizing 
molecular attributes of HN and F as well as the timing, kinetics, 
and interaction between HN and F as a complex in lung tissue 
allows us to consider the HN-F fusion machinery with its various 
components in balance, as a biological complex, tuned for survival 
in the host. 

MATERIALS AND METHODS 

Human airway epithelial cultures. The EpiAirwayAIR-100 system (Mat- 
Tek Corporation) consists of normal, human-derived tracheobronchial 
epithelial cells that have been cultured to form a pseudostratified, highly 
differentiated mucociliary epithelium. Upon receipt from the manufac- 
turer, HAE cultures were handled as previously described (9, 12). 

Monolayer cells. CV- 1 (African green monkey kidney) and 293T (hu- 
man kidney epithelial) cells were grown in Dulbecco's modified Eagle's 
medium (Cellgro; Mediatech) supplemented with 10% fetal bovine serum 
and antibiotics at 37°C in 5% C0 2 . 

Viruses. Titers of HPIV3 virus stocks were assessed by a plaque assay 
(9, 12). Supernatant fluids collected from HAE were used to infect CV-1 
cells. Isolated plaques were enclosed in a 1-mm-diameter plastic cylinder, 
and trypsin was added to disrupt the plaque. Complete medium (100 ul) 
was then added, and the mixture was collected. This fluid was then used to 
infect HAE recurrently to ensure that the plaque retained its morphology 
over 2 to 3 passages for preparation of viral stocks. 

RNA isolation and sequence analysis of HN and F genes. Total RNA 
was isolated from CV1 cells infected with plaque-purified virus or super- 
natant fluids collected from HAE or the HAE cultures themselves, infected 
with 4,000 PFU of HPIV3 WT or the variant viruses. RNA isolation and 
sequencing were performed as previously described, and sequence analy- 
sis was carried out on the DNA products of three separate isolations and 
performed twice on all samples to confirm sequence alterations (24, 25). 

qRT-PCR. To measure the amount of viral genome released, 30 /ul of 
fluid containing released virus was processed for RNA extraction, and 
then the genome was quantified using the real-time one-step qRT-PCR 
pathogen detection kit specific for HPIV3 (PrimerDesign) as per the man- 
ufacturer's instructions. 

Modeling of HN and F and molecular dynamics simulation. For mo- 
lecular dynamics (MD) simulation, the wt or H552Q mutant HN was 
solvated in a water box extending 10 A in all directions beyond the edge of 
the protein, using the TIP3P water model z (26, 27). The system was 
neutralized with a 150 mM concentration of sodium/ chloride ions and 
consisted of a total of 66,000 atoms. Mutagenesis and system preparation 
was performed with the VMD software program (28). MD calculations 
were performed with the NAMD program (29) using the CHARMM27 
force field. The systems were energy minimized with 10,000 steps of con- 
jugate gradient energy minimization, followed by gradual heating from 0 
to 310 K in 30 ps, and then maintained at constant temperature and 
pressure (1.01325 X 10 5 Pa). The simulations were carried out with peri- 
odic box conditions, with a 2-fs time step, a uniform dielectric constant of 
1, a 1 to 4 scaling value of 1 .0, a cutoff of nonbonded forces with a switch- 
ing function starting at 10 A and reaching 0 at 12 A, particle mesh Ewald 
(PME) with a tolerance of 10~ 6 , and all bonds involving hydrogens con- 
strained with the SHAKE algorithm. A production run was performed for 
20 ns, and the trajectories were analyzed with VMD. 

The structure of the prefusion F protein from HPIV3 was homology 
modeled using the HHpred server (http://toollcit.tuebingen.mpg.de 
/hhpred#) and the parainfluenza virus 5 (PIV5) structure as the template 
(18). 

Measuring plaque size. Plaque diameter was measured at a magnifi- 
cation of 1 5 X under a zoom stereomicroscope equipped with a microme- 
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ter. Surface area calculations from 20 plaques were used to establish mean 
plaque areas with standard deviations. 

HN and F constructs and transient expression. Cloned HPIV3 HN 
variant cDNAs were digested with EcoRI and BamHI and ligated into a 
digested pEGFP-C3 vector (BD Biosciences/Clontech). Cloned HPIV3 F 
cDNAs were digested with Xhol and BamHI and ligated into digested 
pCAGGS mammalian expression vectors. Positive clones were sequenced 
to verify the mutations and to ensure that no additional alterations had 
been introduced. Transfections were performed according to the Lipo- 
fectamine Plus or Lipofectamine 2000 manufacturer's protocols (Invitro- 
gen). 

Assay of neuraminidase activity. Monolayers of 293T cells transiently 
expressing WT or mutant HNs were washed once with phosphate- 
buffered saline (PBS) and then incubated for 10 min in 5 mM EDTA in 1 X 
PBS put in pH 5.0 C0 2 -independent medium (Gibco). Cells (150,000) 
were transferred to 96-well plates and incubated at 37°C for 10 min. 2'- 
(4-Methylumbehiferyl)-a-d-N-acetylneuraminic acid sodium salt (MUN 
ANA) (20 mM; Toronto Research Chemical) was added at 1:1 (vol/vol) 
for a final concentration of 10 mM MUNANA. A kinetic reading at 37°C 
was done using a SpectraMax M5 enzyme-linked immunosorbent assay 
(ELISA) reader every 2 min for 1 h. 

Use of receptor-depleted RBCs to assess HN receptor-binding avid- 
ity. Human RBCs were obtained from healthy donors as previously de- 
scribed (24). Partial receptor depletion of RBCs was carried out, and RBCs 
which were partially depleted of their surface sialic acid receptors were 
used to determine the relative receptor binding avidities of variant HN 
molecules as described previously (7, 13). 

Assays for F activation, receptor retention, and receptor release. 
Monolayers of 293T cells transiently expressing F wt or mutant or HN wt 
or mutant were washed three times and incubated with 1 % RBC suspen- 
sions at pH 7.3 for 30 min at 4°C. After rinsing to remove unbound RBCs, 
the plates were transferred to 4, 15, 22, 30, or 37°C for 1 h. After this 
incubation, RBCs that were released to the medium, remained bound, or 
had undergone membrane fusion were quantified as described previously 
(7,13). 

|3-Galactosidase complementation-based fusion assay. The 

j3-galactosidase complementation-based fusion assay was performed as 
previously described (7, 13, 32). 

Cell surface expression assay. Monolayers of 293T cells were trans- 
fected for transient expression of HN or F constructs. Cells were washed 
twice in PBS and then incubated with an anti-HPIV3 HN monoclonal 
antibody supplied by Judy Beeler in 3% bovine serum albumin (BSA), 
0.1% sodium azide in PBS for 1 h. Samples were then washed twice in PBS 
and incubated with 1:100 anti-mouse IgG (H+L) R-phycoerythrin con- 
jugate (Caltag Laboratories). To quantify the amounts of cell surface pro- 
teins in each sample, an indirect immunofluorescence assay was used, 
where beads from the Qifikit kit (Dako) were stained simultaneously with 
samples according to manufacturer's instructions, with the same second- 
ary antibody. Cell surface expression of all samples and beads was per- 
formed with fluorescence-activated cell sorting (FACS) (FACSCalibur; 
Becton, Dickinson) (10, 30, 31). 
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